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Abstract 

Now that the Higgs boson has been observed by the ATLAS and CMS exper- 
iments at the LHC, the next important step would be to measure accurately its 
properties to establish the details of the electroweak symmetry breaking mecha- 
nism. Among the measurements which need to be performed, the determination of 
the Higgs self-coupling in processes where the Higgs boson is produced in pairs is 
of utmost importance. In this paper, we discuss the various processes which allow 
for the measurement of the trilinear Higgs coupling: double Higgs production in 
the gluon fusion, vector boson fusion, double Higgs-strahlung and associated pro- 
duction with a top quark pair. We first evaluate the production cross sections for 
these processes at the LHC with center-of-mass energies ranging from the present 
\/s = 8 TeV to \/s = 100 TeV, and discuss their sensitivity to the trilinear Higgs 
coupling. We include the various higher order QCD radiative corrections, at next- 
to-leading order for gluon and vector boson fusion and at next -to-next-to-leading 
order for associated double Higgs production with a gauge boson. The theoretical 
uncertainties on these cross sections are estimated. Finally, we discuss the various 
channels which could allow for the detection of the double Higgs production signal 
at the LHC and the accuracy on the self-coupling that could be ultimately achieved. 
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1 Introduction 



A bosonic particle with a mass of about 125 GeV has been observed by the ATLAS and 
CMS Collaborations at the LHC [T] and it has, grosso modo, the properties of the long 
sought Higgs particle predicted in the Standard Model (SM) [2]. This closes the first 
chapter of the probing of the mechanism that triggers the breaking of the electroweak 
symmetry and generates the fundamental particle masses. Another, equally important 
chapter is now opening: the precise determination of the properties of the produced 
particle. This is of extreme importance in order to establish that this particle is indeed the 
relic of the mechanism responsible for the electroweak symmetry breaking and, eventually, 
to pin down effects of new physics if additional ingredients beyond those of the SM are 
involved in the symmetry breaking mechanism. To do so, besides measuring the mass, 
the total decay width and the spin-parity quantum numbers of the particle, a precise 
determination of its couplings to fermions and gauge bosons is needed in order to verify 
the fundamental prediction that they are indeed proportional to the particle masses. 
Furthermore, it is necessary to measure the Higgs self-interactions. This is the only way 
to reconstruct the scalar potential of the Higgs doublet field $, that is responsible for 
spontaneous electroweak symmetry breaking, 

Vh = + ^A($t$)2 ; A = ^ and /i^ = -^M^ , (1) 

2 f ^ 2 

with V = 246 GeV. Rewriting the Higgs potential in terms of a physical Higgs boson leads 
to the trilinear Higgs self-coupling Xhhh, which in the SM is uniquely related to the mass 
of the Higgs boson, 

Xhhh = — ■ (2) 

V 

This coupling is only accessible in double Higgs production [3H6]. One thus needs to 
consider the usual channels in which the Higgs boson is produced singly [7], but allows 
for the state to be off mass-shell and to split up into two real Higgs bosons. At hadron 
colliders, four main classes of processes have been advocated for Higgs pair production: 

a) the gluon fusion mechanism, gg — HH, which is mediated by loops of heavy quarks 
(mainly top quarks) that couple strongly to the Higgs boson [5HTT]: 

b) the WW/ZZ fusion processes (VBF), qq' -> V*V*qq' HHqq' {V = W,Z), which 
lead to two Higgs particles and two jets in the final state [5|[T^[T3]: 

c) the double Higgs-strahlung process, qq' V* ^ VHH {V = W, Z), in which the 
Higgs bosons are radiated from either a. W or a. Z boson |14] : 

d) associated production of two Higgs bosons with a top quark pair, pp — )■ tiHH [TB] . 

As they are of higher order in the electroweak coupling and the phase space is small 
due to the production of two heavy particles in the final state, these processes have much 
lower production cross sections, at least two orders of magnitude smaller, compared to the 
single Higgs production case. In addition, besides the diagrams with H* — j- HH splitting, 
there are other topologies which do not involve the trilinear Higgs coupling, e.g. with 
both Higgs bosons radiated from the gauge boson or fermion lines, and which lead to the 
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same final state. These topologies will thus dilute the dependence of the production cross 
sections for double Higgs production on the Xhhh coupling. The measurement of the 
trilinear Higgs coupling is therefore an extremely challenging task and very high collider 
luminosities as well as high energies are required. We should note that to probe the 
quadrilinear Higgs coupling, Xhhhh = 3M|f/t>^, which is further suppressed by a power 
of V compared to the triple Higgs coupling, one needs to consider triple Higgs production 
processes [3l[T6]. As their cross sections are too small to be measurable, these processes 
are not viable in a foreseen future so that the determination of this last coupling seems 
hopeless. 

In Ref. [5], the cross sections for the double Higgs production processes and the 
prospects of extracting the Higgs self-coupling have been discussed for the LHC with 
a 14 TeV center-of-mass (cm.) energy in both the SM and its minimal supersymmetric 
extension (MSSM) where additional channels occur in the various processes. 

In the present paper, we update the previous analysis. In a sense, the task is made 
easier now that the Higgs boson mass is known and can be fixed to Mh ~ 125 GeV. 
However, lower cm. energies have to be considered such as the current one, ^/s = 8 TeV. 
In addition, there are plans to upgrade the LHC which could allow to reach cm. energies 
of about 30 TeV [17] and even up to 100 TeV. These very high energies will be of crucial 
help to probe these processes. 

Another major update compared to Ref. [5] is that we will consider all main processes 
beyond leading order (LO) in perturbation theory, i.e. we will implement the important 
higher order QCD corrections. In the case of the gluon fusion mechanism, gg — > HH, the 
QCD corrections have been calculated at next-to-leading-order (NLO) in the low energy 
limit in Ref. [18]. They turn out to be quite large, almost doubling the production cross 
section at ^/s = 14 TeV, in much the same manner as for single Higgs production [T9] . 
In fact, the QCD corrections for single and double Higgs production are intimately re- 
lated and one should expect, as in the case of gg ^ H, a further increase of the total 
cross section by ~ 30% once the next-to-next-to-leading (NNLO) corrections are also 
included [20l|21]. 

It is well known that for single Higgs production in the vector boson fusion process 
qq' — Hqq' there is no gluon exchange between the two incoming/outgoing quarks as the 
initial and final quarks are in color singlet states at LO. Then the NLO QCD corrections 
consist simply of the known corrections to the structure functions [22142^ . The same can 
be said in the case of double Higgs production qq' — HHqq', and in this paper we will 
implement the NLO QCD corrections to this process in the structure function approach. 
The NNLO corrections in this approach turn out to be negligibly small for single Higgs 
production [25] and we will thus ignore them for double Higgs production. 

In the single Higgs-strahlung process, qq ^ V* ^ VH, the NLO QCD corrections can 
be inferred from those of the Drell-Yan process qq — > V* [26]. This can be extended to 
NNLO [20 | l27 t [28] but, in the case of ZH production, one needs to include the gg initiated 
contribution, gg — )■ ZH [281 [29] cis well as some additional subleading corrections |30] . 
The same is also true for double Higgs-strahlung and we will include in this paper the 
Drell-Yan part of the corrections up to NNLO. In the case of ZHH final states, we will 
determine the additional contribution of the pentagon diagram gg — ZHH which turns 
out to be quite substantial, increasing the total cross section by up to 30% at y/s = 14 
TeV. 
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In the case of the pp — )■ ttHH process, the determination of the cross section at LO is 
aheady rather comphcated. We will not consider any correction beyond this order (that, 
in any case, has not been calculated) and just display the total cross section without 
further analysis. We simply note that the QCD corrections in the single Higgs case, 
pp — > ttH, turn out to be quite modest. At NLO, they are small at ^/s = 8 TeV and 
increase the cross section by less than ^ 20% at y/s = 14 TeV [21]. We also note that 
this channel is plagued by huge QCD backgrounds. 

In addition, the electroweak radiative corrections to these double Higgs production 
processes have not been calculated yet. Nevertheless, one expects that they are similar 
in size to those affecting the single Higgs production case, which are at the few percent 
level at the presently planned LHC cm. energies [321133] • They should thus not affect the 
cross sections in a significant way and we will ignore this issue in our analysis. 

After determining the i^-factors, i.e. the ratios of the higher order to the lowest 
order cross sections consistently evaluated with the value of the strong coupling as and 
the parton distribution functions taken at the considered perturbative order, a next step 
will be to estimate the theoretical uncertainties on the production cross sections in the 
various processes. These stem from the variation of the renormalization and factorization 
scales that enter the processes (and which gives a rough measure of the missing higher 
order contributions), the uncertainties in the parton distribution functions (PDF) and 
the associated one on the strong coupling constant and, in the case of the gg — HH 
process, the uncertainty from the use of an effective approach with an infinitely heavy 
virtual top quark, to derive the NLO corrections (see also Ref. [31] )• This will be done in 
much the same way as for the more widely studied single Higgs production case [33] 135]. 

Finally, we perform a preliminary analysis of the various channels in which the Higgs 
pair can be observed at the LHC with a cm. energy of y/s = 14 TeV assuming up 
to 3000 fb"^ collected data, and explore their potential to probe the Xhhh coupling. 
Restricting ouselves to the dominant gg — ?■ HH mechanism in a parton level approachu, 
we first discuss the kinematics of the process, in particular the transverse momentum 
distribution of the Higgs bosons and their rapidity distribution at leading order. We then 
evaluate the possible cross sections for both the signal and the major backgrounds. As the 
Higgs boson of a mass ~ 125 GeV dominantly decays into 6-quark pairs with a branching 
ratio of ~ 60% and other decay modes such as if — )■ 77 and H — > WW* — 2i2i> are 
rare [H], and as the production cross sections are already low, we will focus on the three 
possibly promising detection channels gg — t- HH — )■ 6677, bbrT and bbW~^W~. Very high 
luminosities, C(ab~^) would be required to have some sensitivity on the Xhhh coupling. 

The rest of the paper is organized as follows. In the next section, we discuss the 
QCD radiative corrections to double Higgs production in the gluon fusion, vector boson 
fusion and Higgs-strahlung processes (the tiHH process will be only considered at tree- 
level) and how they are implemented in the programs HPAIR [12], VBFNLO [13] and a code 
developed by us to evaluate the inclusive cross sections in Higgs-strahlung processes. In 

^ Early and more recent parton level analyses of various detection channels have been performed in 
Refs. [36ll38] with the recent ones heavily relying on jet-substructure techniques [39]. However, a full 
and realistic assessment of the LHC to probe the trilinear coupling would require the knowledge of the 
exact experimental conditions with very high luminosities and a full simulation of the detectors which is 
beyond the scope of this paper. Only the ATLAS and CMS Collaborations are in a position to perform 
such detailed investigations and preliminary studies have already started |40j . 
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section 3, we evaluate the various theoretical uncertainties affecting these cross sections 
and collect at Mh = 125 GeV the double Higgs production cross sections at the various 
LHC energies. We also study the sensitivity in the different processes to the trilinear 
Higgs self-coupling. Section 4 will be devoted to a general discussion of the channels that 
could allow for the detection of the two Higgs boson signal at a high-luminosity 14 TeV 
LHC, concentrating on the dominant gg — ?■ HH process, together with an analysis of the 
major backgrounds. A short conclusion is given in a last section. 



2 Higgs pairs at higher orders in QCD 

Generic Feynman diagrams for the four main classes of processes leading to double Higgs 
production at hadron colliders, gluon fusion, WW/ZZ fusion, double Higgs-strahlung 
and associated production with a top quark pair, are shown in Fig. 1. As can be seen in 
each process, one of the Feynman diagrams involves the trilinear Higgs boson coupling, 
Xhhh = 3M^/w, which can thus be probed in principle. The other diagrams involve 
the couplings of the Higgs boson to fermions and gauge bosons and are probed in the 
processes where the Higgs particle is produced singly. 

In this section we will discuss the production cross sections for the first three classes 
of processes, including the higher order QCD corrections. We will first review the gluon 
channel and then we will move on to the higher-order corrections in the weak boson fusion 
and Higgs-strahlung channels. 



2.1 The gluon fusion process 

The gluon fusion process is - in analogy to single Higgs production - the dominant Higgs 
pair production process. The cross section is about one order of magnitude larger than 
the second largest process which is vector boson fusion. As can be inferred from Fig. [1^) 
it is mediated by loops of heavy quarks which in the SM are mainly top quarks. Bottom 
quark loops contribute to the total cross section with less than 1% at LO. 

The process is known at NLO QCD in an effective field theory (EFT) approximation 
by applying the low energy theorem (LET) [T ^ H H H5] which means that effective cou- 
plings of the gluons to the Higgs bosons are obtained by using the infinite quark mass 
approximation. The hadronic cross section at LO is given by 



f (It (Tlo(s = Ts) j —fg{x;f-il)fg 



(3) 



with s being the hadronic cm. energy, Tq = ^Mj^/s, and fg the gluon distribution function 
taken at a typical scale fip specified below. The partonic cross section at LO, (Tlo, can 
be cast into the form 



^Loigg -> HH) 



dt 



256(27r)3 



>^HHH V 



Ml 
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(4) 



where 



Ml 
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H 



(5) 
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(a) gg double-Higgs fusion: gg — HH 
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(b) WW/ZZ double-Higgs fusion: qq' — )■ HHqq' 




(c) Double Higgs-strahlung: qq' ZHH/WHH 





(d) Associated production with top-quarks: qq/ gg — > ttHH 
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Figure 1: Some generic Feynman diagrams contributing to Higgs pair production at hadron 
colliders. 

with s and t denoting the partonic Mandelstam variables. The triangular and box form 
factors -Fa, Fb and Gu approach constant values in the infinite top quark mass limit, 



3' 



3' 



(6) 



The expressions with the complete mass dependence are rather lenghty and can be found 
in Ref. [H] as well as the NLO QCD corrections in the LET approximation in Ref. [18] . 

The full LO expressions for -FA,-Pb and Gu are used wherever they appear in the 
NLO corrections in order to improve the perturbative results, similar to what has been 
done in the single Higgs production case where using the exact LO expression reduces the 
disagreement between the full NLO result and the LET result [TtfTU]. 

For the numerical evaluation we have used the publically available code HPAIR [32] in 
which the known NLO corrections are implemented. As a central scale for this process 
we choose 

AtO = = jJ'F = Mhh , (7) 

where Mhh denotes the invariant mass of the Higgs boson pair. The i^-factor, describing 
the ratio of the cross section at NLO using NLO PDFs and NLO as to the leading order 
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cross section consistently evaluated with LO PDFs and LO ag, for this process is 

K ~ 2.0 (1.5) for = 8 (100) TeV . (8) 

2.2 The vector boson fusion process 

The structure of the Higgs pair production process through vector boson fusion [12] is 
very similar to the single Higgs production case. The vector boson fusion process can 
be viewed as the double elastic scattering of two (anti) quarks with two Higgs bosons 
radiated off the weak bosons that fuse. In particular this means that the interference 
with the double Higgs-strahlung process qq' — >■ V*HH — qq'HH is negligible and this 
latter process is treated separately. This is justified by the kinematics of the process 
with two widely separated quark jets of high invariant mass and by the color flow of 
the process. This leads to the structure function approach that has been applied with 
success to calculate the QCD corrections in the vector boson fusion production of a single 
Higgs boson [22] - [25] . The diagrams contributing at NLO QCD order are shown in Fig. [2l 
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Figure 2: Generic diagrams contributing to the NLO corrections to qq' HHqq' . Shown 
are the LO diagram (upper left) and the NLO corrections for the upper quark line. The 
blob of the VVHH vertex is a shortcut for the three diagrams shown in Figs. [Tl)) and\^ 

For simplicity only the diagrams with the QCD corrections to the upper quark line are 
shown. The calculation involving the second quark line is identical. The blob of the 
vertex VVHH is a shortcut for the diagrams depicted in Fig. [3], which include charged 
currents (CC) with W"^ bosons and neutral currents (NC) with a Z boson exchange. As 
can be seen only one of the three diagrams involves the trilinear Higgs coupling. The 
other diagrams act as irreducible background and lower the sensitivity of the production 
process to the Higgs self-coupling. 

We have calculated the NLO QCD corrections in complete analogy to the single Higgs 
VBF process [23] • The real emission contributions are given by a gluon attached to the 
quark lines either in the inital or the final state and from the gluon-quark initial state. As 
we are working in the structure function approach, the corrections of the upper and lower 
quark lines do not interfere and are simply added incoherently. The amplitudes have the 
following structure, 

AHH,,'^n':v*JiJi . (9) 
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Figure 3: Diagrams contributing to the VVHH vertex. 



where Tyiy, stands for the tensor structure of the diagrams depicted in Fig.[3]and J^'^' are 
the quark currents of the upper and lower hues, respectively, with four-momenta g, g'. The 
calculation is done numerically using the Catani- Seymour dipole subtraction method [16] 
to regularize the infrared divergencies. The formulae for the subtraction terms as well as 
the finite corrections are identical to the ones for single Higgs VBF production as only 
the quark currents are involved. They can be found in Ref. |23j . 

We have implemented this calculation in the VBFNLO code [13] in which we have 
provided the tensor structure depicted in Fig. [3] which has been calculated with Mad- 
Graph [1?]. Up to now the VBFNLO implementation only involves on-shell Higgs pairs. 
We have found an increase of ~ +7% of the total cross section compared to the LO result 
when using the central scale 

fJ'O = fJ'R = fJ'F = Qv*, (10) 

with Qv* being the momentum of the exchanged weak bosons This 
result is in agreement with a previous calculation done in the context of the two Higgs 
doublet model 



2.3 The Higgs— strahlung process 

The production of a Higgs pair in association with a vector boson has been calculated for 
the first time quite a while ago [H] and shares common aspects with the single Higgs- 
strahlung process. The NLO corrections can be implemented in complete analogy to 
single Higgs-strahlung [2H]- We will update in this paper the former results and present 
the NNLO corrections to the WHH and ZHH inclusive production cross sections. These 
calculations have been implemented in a code which shall become publically available. 

At LO the process pp — VHH {V = W, Z) is given by quark-antiquark annihilations 
in s-channel mediated processes involving three Feynman diagrams, see Fig. [Tt). As 
can be seen only one of the three diagrams involves the trilinear Higgs coupling. The 
sensitivity to this coupling is then diluted by the remaining diagrams. After integrating 
over the azimuthal angle we are left with the following partonic differential cross section 
with s being the partonic cm. energy (see also Ref. [5]), 



dxidx2 1U9V2tt^s{1 - 

h 



+ 



1 



8 4/iy(l - a;i + /iH - /iyj 



X 



+ 



/2 



1 — Xi + flH — f^V ^ — X2 + Hh — f^V 



^In order to stay within the perturbative regime a cut Qy > 2 GeV has to be imposed, see Ref. [22] . 
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where we use of the following notation, 



2Eh 2Ev 



and the reduced couplings of the quarks to the vector bosons, ag = Vq = a/2 for = ly 
and any quark g, a„ = 1 and f „ = 1 — 8/3 sin^ 9^ for q = u,s and V = Z , aa = —1 and 
Vd = —1 + 4/3 sin^ 9w for q = d,c,b and V = Z. The coefficients /j as well as Chhh are 

/o = ^J'V [(2 - xi - X2)^ + 8/iv/] , 

/i = (/iy - 1 + xi)^ - (1 - a;i) (1 - xi + /iy - /iyxi - 4/iy) 

(/iy - 4;Uh) (1 - 4/1^) - /U^, 
/2 = (2/iy + Xi + X2) [^^V (Xi + X2 - 1 + /iy - 8/iH) 

- (1 - xi) (1 - X2) (1 + ^lv)] + (1 - xi)' (1 - X2)' 

+ (1 - Xi) (1 - X2) [/iy + 1 + ^^lH (1 + ^^v)\ 
+A^,H^iv (1 + /iy + 4/i//) + /iy, 

/s = Xi (Xi - 1) (//y + Xi - 1) - (1 - X2) (2 - Xi) (1 - Xi + /iy) 
+2/iy (/iy + 1 - 4/i^) , 

^ ''^ ^HHH . 2 



^HHH - TJ2 



My a;i + X2 - 1 + /iy - //// I - xi + jiH - ^-v 

1 — X2 + /i/f — /iy /iy 

The coefficient Chhh includes the trilinear Higgs coupling Xhhh- 

In order to obtain the full hadronic section, the differential partonic cross section of 
Eq. ( ITT]) is convoluted with the quark parton distribution functions, fqt taken at a 
typical scale ixp specified below: 

a{pp l^M) = XI / / —fi(^'^ ^^F) fq/q' (-; A^f) ^vhh{s = Ts) , (14) 



q,q' 



where s stands for the hadronic cm. energy and tq = {2Mh + My)"^ / s. The total partonic 
cross section (Tvhh has been obtained after the integration of Eq. lITT]) over Xi,X2. 

The calculation of the NLO QCD corrections is similar to the single Higgs-strahlung 
case. In fact, this process can be viewed as the Drell-Yan production pp — V* followed 
by the splitting process V* — ?■ VHH. The off-shell vector boson can have any momentum 
k"^ with {2Mh + My)^ < < s. This factorization is in principle valid at all orders for 
the Drell-Yan like contributions and leads, after folding with the PDF, to 

/"i r/Z'*-' 

a{pp ^ VHH) = drj^-^ dza^^{zTs)Aij{ij -> V*) , (15) 



with 

dr 



|V.(a:;/i|)/,(^;/i|)^. (16) 
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In the expressions above ij stands for any initial partonic subprocess, Aj^ is the Drell- 
Yan correction, z = k'^/s and o"^*^ stands for the LO partonic cross section of the process 
qq' ->V*-^ VHH. At LO we have = 5iq5jqiqi5{l — z). In Fig.lHthe generic diagrams 
which contribute at NLO to the Drell-Yan process qq' — )■ V* are depicted. The NLO QCD 
corrections increase the total cross section by ~ +17% at 14 TeV for Mh = 125 GeV. 




Figure 4- Feynman diagrams contributing to the NLO QCD corrections for Drell-Yan 
production. 

We have calculated the NNLO corrections, which have not been available so far, in 
the same way except for the process involving a Z boson. In fact there are additional 
contributions that are specific to the case of a Z boson, involving an effective Zgg vertex. 
Similar to what is stated in Ref. [2S] for the single Higgs production case, only the specific 
gluon fusion initiated process will be of non-negligible contribution and will be described 
below. Let us start with the NNLO QCD Drell-Yan contribution. Some generic diagrams 
contributing to the NNLO corrections to qq' — V* are shown in Fig. [5l We apply the 
procedure as described by Eq. f[T^ and the expression is then given by 




Figure 5: Some Feynman diagrams contributing at NNLO QCD to Drell-Yan production. 
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with 5vz = 1(0) for V = Z{W) and 




The expressions for the coefficients A^'^^^''^\z) refer to the NLO and NNLO corrections, 
respectively. As they are too lengthy to be reproduced here, we refer the reader to the 
appendix B of Ref. [27] and to Ref. [20]. The expressions given there have to be rescaled 
by a factor of (vr/a^)*, and M = fip, R = hr. In our calculation we have included the 
full CKM matrix elements in the quark luminosity as well as the initial bottom quark 
contribution. We use the central scale 

fJ'O = m = fJ'F = MvHH , (19) 

where Mvhh denotes the invariant mass of the VHH system. 

The Drell-Yan NNLO QCD corrections Eq. ( |T8l) turn out to be very small. They 
typically increase the cross section by a few percent at 14 TeV. 

The last contribution Aagg^zHH, see diagrams in Fig. [6l is only present in the case 
of Higgs pair production in association with a Z boson. It stems from the process 
gg — 7- ZHH which is loop-mediated already at LO. Being of order it contributes 
to the total cross section pp — )■ ZHH at NNLO QCD. The process is mediated by quark 
loops. Furthermore, there are also box and pentagon topologies where again only top and 
bottom quarks contribute as the Yukawa couplings to light quarks can be neglected. At 
the LHC the contribution of the gluon fusion channel is substantial in contrast to the sin- 
gle Higgs production case. Indeed, while in the latter the contribution is of order ~ +5% 
compared to the NNLO QCD Drell-Yan contribution, in the case of Higgs pair production 
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Figure 6: Some generic diagrams contributing to gg — j- ZHH. For the triangle+hox 
topologies, only those involving the trilinear Higgs couplings are depicted. 



it contributes with ~ +20 ■ ■ - + 30% depending on the cm. energy. This enhancement can 
be explained by the additional pentagon topology which a) involves two top Yukawa cou- 
plings and b) softens the destructive interference between the triangle and box diagrams 
that is present in the single Higgs production case. Furthermore, the suppression by a 
power (as/vr)^ is partly compensated by the increased gluon luminosity at high energies. 
This explains why this channel, which has been calculated using FeynArts/FormCalc |49] . 
should be taken into account. It also implies that the scale variation in pp — ZHH will 
be worse than in pp — )■ WHH because of the LO gluon fusion mechanism appearing at 
NNLO. 



3 Cross Sections and sensitivity at the LHC 

In this section we will present the results for the calculation of the total cross sections 
including the higher-order corrections discussed in the previous section as well as the 
various related theoretical uncertainties. We will use the MSTW2008 PDF set [SD] as our 
reference set. We choose the following values for the W , Z and top quark masses and for 
the strong coupling constant at LO, NLO and NNLO, 

Mh/ = 80.398 GeV, = 9L1876 GeV, = 173.1 GeV, 

aLO(M|) = 0.13939, afLO(M|) = 0.12018, a^NLO(M|) = 0.11707. 

The electromagnetic constant a is calculated in the scheme from the values of Mw 
and Mz given above. For the estimate of the residual theoretical uncertainties in the 
various Higgs pair production processes we considered the following uncertainties: 

1. the scale uncertainty, stemming from the missing higher order contributions and 
estimated by varying the renormalization scale /ir and the factorization scale 

in the interval |yUo < fiRyfip < 2/io with some restrictions on the ratio (j^r/ fip 
depending on the process; 

2. the PDF and related errors. The PDFs are non-perturbative quantities fitted 
from the data and not calculated from QCD first principles. It is then compulsory 
to estimate the impact of the uncertainties on this fit and on the value of the strong 
coupling constant a<((M|) which is also fitted together with the PDFs; 

3. in the case of the gluon fusion process there is a third source of uncertainties which 
comes from the use of the effective field theory approximation to calculate the QCD 
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NLO corrections, where top loops are taken into account in the infinite top mass 
approximation and bottom loops are neglected. 

In the following we will present results for Mh = 125 GeV. Note that the results for the 
total cross sections and uncertainties are nearly the same for Mh = 126 GeV. The total 
cross sections at the LHC for the four classes of Higgs pair production processes are shown 
in Fig. [7] as a function of the cm. energy. For all processes the numerical uncertainties 
are below the permille level and have been ignored. The central scales which have been 
used are {hr = fip = /io) 




= Mt + MH. (20) 



0.1 I ' ' ' ' 1 
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Figure 7: The total cross sections including higher-order correction at the LHC for Higgs 
pair production in the main channels - gluon fusion (red/full), VBF (green/ dashed), 
Higgs- strahlung (blue/dotted), associated production with ti (violet/ dotted with small dots) 
- as a function of the cm. energy with Mh = 125 GeV. The MSTW2008 PDF set has 
been used and higher-order corrections are included as discussed in section 2. 

As can be inferred from the figure and also seen in Table [1] the largest cross section is 
given by the gluon fusion channel which is one order of magnitude larger than the vector 
boson fusion cross section. All processes are ~ 1000 times smaller than the corresponding 
single Higgs production channels, implying that high luminosities are required to probe 
the Higgs pair production channels at the LHC. 

3.1 Theoretical uncertainties in the gluon channel 

3.1.1 Theoretical uncertainty due to missing higher order corrections 

The large i^-factor for this process of about 1.5 — 2 depending on the cm. energy 
shows that the inclusion of higher order corrections is essential. An estimate on the size 
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x/i [TeV] 




NLO T-fT-il 
'^qq'^HHqq' i^"J 


^qq'^WHH i^"J 




^q?/aa-,ttHH [fb] 


8 


8.16 


0.49 


0.21 


0.14 


0.22 


14 


33.89 


2.01 


0.57 


0.42 


1.09 


33 


207.29 


12.05 


1.99 


1.68 


8.37 


100 


1417.83 


79.55 


8.00 


8.27 


82.69 



Table 1: The total Higgs pair production cross sections in the main channels at the LHC 
(in fb) for given cm. energies (in TeV) with Mh = 125 GeV. The central scales which 
have been used are described in the text. 



of the uncertainties due to the missing higher order corrections can be obtained by a 
variation of the factorization and renormahzation scales of this process. In analogy to 
single Higgs production studies [33ll35] we have estimated the error due to missing higher 
order corrections by varying ^n, fip in the interval 

^/io < At/? = < 2/io. (21) 

As can be seen in Fig.[8]we find sizeable scale uncertainties of order ~ +20%/— 17% 
at 8 TeV down to +12%/ — 10% at 100 TeV. Compared to the single Higgs production 
case the scale uncertainty is twice as large [33l[35] . However, this should not be a surprise 
as there are NNLO QCD corrections available for the top loop (in a heavy top mass 
expansion) in the process gg ^ H while they are unknown for the process gg — )■ HH. 
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Figure 8: Scale uncertainty for a scale variation in the interval |/io < I^r = I^f ^ 2/io in 
^{.99 ~^ HH) at the LHC as function of y/s at Mh = 125 GeV. In the insert the relative 
deviations to the results for the central value fiQ = fiji = fip = Mhh ore shown. 

3.1.2 PDF and as errors 

The parametrization of the parton distribution functions is another source of theoretical 
uncertainty. First there are pure theoretical uncertainties coming from the assumptions 
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Figure 9: The NLO cross section (j{gg — )■ HH + X) at the LHC as a function of the cm. 
energy for = 125 GeV, when using different NLO PDF sets. In the insert the cross 
sections normalized to the cross section calculated with the MSTW PDF set are shown. 



made on the parametrization, e.g. the choice of the parametrization, the set of input 
parameters used, etc. Such uncertainties are rather difficult to quantify. A possibihty 
might be to compare different parameter sets, such as MSTW [50], CTIO [51], ABMll [52], 
GJR08 [S3], HERA 1.5 [M] and NNPDF 2.3 [33]. This is exemphfied in Fig. where 
the predictions using the six previous PDF sets are displayed. As can be seen there 
are large discrepancies over the wole considered cm. energy range. At low energies 
the smallest prediction comes from ABMll which is ~ 22% smaller than the prediction 
made with the MSTW set while at high energies ABMll and MSTW lead to similar 
results whereas the result obtained with the GJR08 set deviates by ~ —15%. The CTIO 
predictions show about —5% difference over the whole energy range with respect to the 
cross section obtained with MSTW while the HERA prediction starts from lower values 
and eventually reaches the CTIO result. Finally the cross sections calculated with the 
NNPDF set decrease over the energy range, starting from being very similar to the MSTW 
result to reach at a/s = 100 TeV the one calculated with CTIO. 

Another source of uncertainty due to the PDF sets comes from the experimental 
uncertainties on the fitted data. The so-called Hessian method, used by the MSTW 
collaboration, provides additional PDF sets next to the best- fit PDF. Additional 2Np£)F 
sets refiect the ±lcr variation around the minimal of all Nppjp parameters that enter 
the fit. Using the 90% CL error PDF sets provided by the MSTW collaboration a PDF 
error of about 6% is obtained for i/s = 8 TeV. The uncertainty shrinks to ~ 2% for 

= 100 TeV. 

In addition to the PDF uncertainties, there is also an uncertainty due to the errors 
on the value of the strong coupling constant a^. The MSTW collaboration provides 
additional PDF sets such that the combined PDE+a^j uncertainties can be evaluated [56] . 
At NLO the MSTW PDF set uses 

as{Mz) = 0.12018+°:°°}22^g^^ gg^ +0.00317 g^a/^ . (22) 

As the LO process is already C(a^), uncertainties in as can be quite substantial. 
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The combined PDF and as error is much larger than the pure PDF error. At 8 TeV 
the PDF error of +5.8%/-6.0% rises to a combined error of +8.5%/-8.3%. At 33 TeV the 
rise is even larger - from the pure PDF error of +2.5%/— 2.7% to the combined PDF+a^ 
error of +6.2%/ -5.4%. 

There is also a theoretical uncertainty on as stemming from scale variation or ambi- 
guities in the heavy flavour scheme definition. The MSTW collaboration estimates this 
uncertainty for as to Aas{Mz) = ±0.003 [56]. However this uncertainty is already in- 
cluded in the scale uncertainty on the input data sets used to fit the PDF and has been 
taken into account in the definition of the PDF+a^ uncertainty. Thus, it does not have 
to be taken into account separately. 

However, even if all these uncertainties for the MSTW PDF set are taken into ac- 
count, the different PDF set predictions do not agree. There might be agreement if also 
uncertainties of the other PDF sets are taken into account, as done in Ref. [35] for the 
case of single Higgs production. 

3.1.3 The effective theory approach 

The last source of theoretical errors that we consider is the use of the LET for the calcu- 
lation of the NLO corrections. At LO it was found in Ref. [1] that for ^/i = 16 TeV the 
LET underestimates the cross section by (9(20%). Furthermore this can be even worse for 
different energies, not to mention the fact that the LET approximation produces incorrect 
kinematic distributions [36] • The reason is that the LET is an expansion in m^p ^ \f§. 
Such an expansion works very well for single Higgs production, since = Mh (at LO) 
for the production of an on-shell Higgs boson whereas in Higgs pair production we have 

2Mh <VI<Vs. (23) 

This means that for Higgs pair production rritop ^ is never fulfilled for Mjj = 125 GeV 
so that the LET approximation is not valid at LO j34j . 

At NLO, however, the LET approximation works much better in case the LO cross 
section includes the full mass dependence. The reason is that the NLO corrections are 
dominated by soft and collinear gluon effects. They factorize in the Born term and in the 
NLO correction contributions, meaning that the /^-factor is not strongly affected from 
any finite mass effects. Based on the results for the single Higgs case [19] we estimate 
the error from applying an effective field theory approach for the calculation of the NLO 
corrections conservatively to 10%. 

3.1.4 Total uncertainty 

In order to obtain the total uncertainty we follow the procedure advocated in Ref. |57] . 
Since quadratic addition is too optimistic (as stated by the LHC Higgs Cross Section 
Working Group, see Ref. [33]), and as the linear uncertainty might be too conservative, the 
procedure adopted is a compromise between these two ways of combining the individual 
theoretical uncertainties. We first calculate the scale uncertainty and then add on top of 
that the PDE+a^ uncertainty calculated for the minimal and maximal cross sections with 
respect to the scale variation. The LET error is eventually added linearly. This is shown 
in Fig. [To] where we display the total cross section including the combined theoretical 
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Figure 10: The total cross section (black/ full) of the process gg HH + X at the 
LHC for Mh = 125 GeV as a function of y/s including the total theoretical uncertainty 
(red/dashed). The insert shows the relative deviation from the central cross section. 



uncertainty. It is found to be sizeable, ranging from ~ +42%/ — 33% at 8 TeV down to 
~ +30%/ — 25% at 100 TeV. The numbers can be found in Table [2l 



[TeV] 




Scale 


[%] 


PDF 


[%] 


PDF+a, 


[%] 


EFT [%] 


Total 


[%] 


8 


8.16 


+20.4 


-16.6 


+5.8 


-6.0 


+8.5 - 


-8.3 


±10.0 


+41.5 


-33.3 


14 


33.89 


+18.2 


-14.7 


+3.9 


-4.0 


+7.0 - 


-6.2 


±10.0 


+37.2 


-29.8 


33 


207.29 


+15.2 


-12.4 


+2.5 


-2.7 


+6.2 - 


-5.4 


±10.0 


+33.0 


-26.7 


100 


1417.83 


+12.2 


-9.9 


+2.0 


-2.7 


+6.2 - 


-5.7 


±10.0 


+29.7 


-24.7 



Table 2: The total Higgs pair production cross section at NLO in the gluon fusion process 
at the LHC (in fb) for given cm. energies (in TeV) at the central scale fJ^F = f^R = Mhh, 
for Mh = 125 Ge V. The corresponding shifts due to the theoretical uncertainties from the 
various sources discussed are shown as well as the total uncertainty when all errors are 
added as described in the text. 



3.2 VBF and Higgs— strahlung processes 

We will not repeat the detailed description of the previous uncertainties in this subsection 
and only summarize how they affect the VBF and Higgs-strahlung inclusive cross sections. 
In both channels, only the scale uncertainties and the PDF+a^ errors are taken into 
account, the calculation being exact at a given order. 

3.2.1 The VBF channel 

As stated in section 2 we use the central scale fio = fJ^R = fJ^F = Qv* ? that is the momentum 
transfer of the exchanged weak boson. Note that a cut of Qv* > 2 GeV has to be applied 
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as stated in the previous section. The scale uncertainty is calculated in exactly the same 
way as for the gluon fusion mechanism, exploring the range /io/2 < fiR,fiF < 2/io. We 
have checked that imposing the restriction 1/2 < fin/fip < 2 does not modify the final 
result. We obtain very small scale uncertainties ranging from ~ ±2% at 8 TeV down to 
~ ±1% and even lower at 33 TeV as can be seen in Fig. [11] (left). This is in sharp contrast 
with the ±8% uncertainty obtained at LO at 8 TeV for example, which illustrates the 
high level of precision already obtained with NLO QCD corrections. 

The 90% CL PDF+a^ uncertainties are calculated following the recipe presented in the 
gluon fusion subsection. The PDF+a^ dominates the total error, ranging from +7%/— 4% 
at 8 TeV down to +5%/ -4% at 33 TeV. 
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Figure 11: Scale uncertainty for a scale variation in the interval |/io < fiR, fip < 2/io 
(left) and total uncertainty bands (right) in a{qq' — t- HHqq') at the LHC as a function 
of at Mfj = 125 GeV. The inserts show the relative deviations to the cross section 
evaluated at the central value fio = fiR = fip = Qv* ■ 



The total error has been obtained by adding linearly the scale and PDF+Os uncer- 
tainties, given the small variation of the cross section with respect to the choice of the 
scale. This process has a total theoretical uncertainty which is always below 10%, from 
+9%/-6% at 8 TeV to +6%/-4% both at 33 and 100 TeV as can be read off Table [3] 
The total uncertainty is displayed in Fig. [TT] (right) as a function of the cm. energy. The 
QCD corrections drastically reduce the residual theoretical uncertainty. 



3.2.2 The associated Higgs pair production with a vector boson 

The cross section is calculated with the central scale fio = fJ^R = fip = Mvhh which is 
the invariant mass of the W/Z + Higgs pair system. The scales are varied in the interval 
Aio/2 < fJ'R = fJ'F < 2yUo- The factorization and renormalization scales can be chosen 
to be the same as the impact of taking them differently is totally negligible, given the 
fact that the scale fiR only appears from NLO on and that we have a NNLO calculation 
which then reduces any non-negligible contribution arising from the difference between 
renormalization and factorization scales. As noticed in section 2, the scale uncertainty is 
expected to be worse in the ZHH channel because of the significant impact of the gluon 
fusion contribution. This is indeed the case as we obtain a scale uncertainty below ±0.5% 
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^/i [TeV] 


[fb] 


Scale 


[%] 


PDF 


[%] 


PDF+a, [%] 


Total 


[%] 


8 


0.49 


+2.3 


-2.0 


+5.2 


-4.4 


+6.7 -4.4 


+9.0 


-6.4 


14 


2.01 


+1.7 


-1.1 


+4.6 


-4.1 


+5.9 -4.1 


+7.6 


-5.1 


33 


12.05 


+0.9 


-0.5 


+4.0 


-3.7 


+5.2 -3.7 


+6.1 


-4.2 


100 


79.55 


+1.0 


-0.9 


+3.5 


-3.2 


+5.2 -3.2 


+6.2 


-4.1 



Table 3: The total Higgs pair production cross section at NLO in the vector boson fusion 
process at the LHC (in fb) for given cm. energies (in TeV) at a central scale fip = f^R = 
Qv* for Mh = 125 GeV. The corresponding shifts due to the theoretical uncertainties 
from the various sources discussed are also shown as well as the total uncertainty when 
all errors are added linearly. 



in the WHH channel whereas the uncertainty in the ZHH channel is A'^ ~ ±3% over 
the entire energy range (slighty more at higher energies) as can be seen in Fig. [T2] (left). 

The total PDF+ag error that we obtain is very similar for the two channels WHH 
and ZHH. It varies from ~ ±4% at 8 TeV down to ~ ±3% at 33 TeV, with a slightly 
higher uncertainty at 100 TeV. 

The total error has been obtained exactly as in the VBF case, given the very small 
variation of the cross section with respect to the scale choice. It is dominated by the 
PDF+a, uncertainty and amounts to +5%/^% (+4%/^%) at 8 (100) TeV in the WHH 
channel and +7%/ -5% (+8%/ -8%) at 8 (100) TeV in the ZHH channel. The total 
theoretical uncertainty in the Higgs-strahlung channels is less than 10%. The numbers are 
given in Tables H] and |5l The total uncertainty bands for the WHH and ZHH channels 
are displayed in Fig. [T2] (right). 
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Figure 12: Scale uncertainty for a scale variation in the interval |/io < /^r, < 2/io (left) 
and total uncertainty bands (right) in Higgs pair production through Higgs-strahlung at 
NNLO QCD at the LHC as a function of ^/s for Mh = 125 CeV. The inserts show the 
relative deviations to the cross section evaluated with the central scale fJ^o = fJ^R = fip = 

MvHH- 
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^/i [TeV] 




Scale 


[%] 


PDF 


[%] 


PDF+a, [%] 


Total 


[%] 


8 


0.21 


+0.4 


-0.5 


+4.3 


-3.4 


+4.3 -3.4 


+4.7 


-4.0 


14 


0.57 


+0.1 


-0.3 


+3.6 


-2.9 


+3.6 -3.0 


+3.7 


-3.3 


33 


1.99 


+0.1 


-0.1 


+2.9 


-2.5 


+3.4 -3.0 


+3.5 


-3.1 


100 


8.00 


+0.3 


-0.3 


+2.7 


-2.7 


+3.8 -3.4 


+4.2 


-3.7 



Table 4- The total Higgs pair production cross sections at NNLO in the qq' — WHH 
process at the LHC (in fb) for different cm. energies (in TeV) at the central scale 
IJ^F = iJ^R = My/HH for Mfj = 125 GeV. The corresponding shifts due to the theoreti- 
cal uncertainties from the various sources discussed are also shown as well as the total 
uncertainty when all errors are added linearly. 



V7s [TeV] 


-ThW m 


Scale 


[%] 


PDF 


[%] 


PDF+a, 


[%] 


Total 


[%] 


8 


0.14 


+3.0 


-2.2 


+3.8 


-3.0 


+3.8 - 


-3.0 


+6.8 


-5.3 


14 


0.42 


+4.0 


-2.9 


+2.8 


-2.3 


+3.0 - 


-2.6 


+7.0 


-5.5 


33 


1.68 


+5.1 


-4.1 


+1.9 


-1.5 


+2.7 - 


-2.6 


+7.9 


-6.7 


100 


8.27 


+5.2 


-4.7 


+1.9 


-2.1 


+3.2 - 


-3.2 


+8.4 


-8.0 



Table 5: Same as Table\^for ZHH production using the central scale fJ^F = fJ'R = Mzhh- 

3.3 Sensitivity to the trilinear Higgs coupling in the main chan- 
nels 

We end this section by a brief study of the sensitivity of the three main channels to 
the trihnear Higgs couphng that we want to probe. Indeed, as can be seen in Fig. [H 
all processes do not only involve a diagram with the trilinear Higgs couplings but also 
additional contributions which then dilute the sensitivity. In order to study the sensitivity 
within the SM, we rescale the coupling Xhhh in terms of the SM trilinear Higgs coupling 
as Xhhh = i^X^fjj^j. 

In Fig. [13] the three main Higgs pair production cross sections are displayed as a 
function of n for three different cm. energies ^/s = 8, 14 and 33 TeV. The left panels 
show the total cross sections while the right panels show the ratio between the cross 
sections at a given Xhhh = f^X^^ and the SM cross section with k = 1. 

As can be seen, the most sensitive channel is by far the VBF production mode, in 
particular for low and high values of k. The shapes of the cross sections with respect to 
a variation of k are the same in all channels and at all energies with a minimum reached 
at K ~ — 1, 2 and 3 for Higgs-strahlung, VBF and gluon fusion, respectively. The right 
panels of Fig. [13] also show that the sensitivity decreases when a/s increases. This is to 
be expected as the diagrams involving the trilinear Higgs self-coupling are mediated by 
s-channel propagators which get suppressed with increasing energy, so that the relative 
importance of these diagrams with respect to the remaining ones is suppressed. 

It is important to note that the sensitivity tested here does not give information on the 
sensitivity to Higgs self-couplings in models beyond the SM. It only tests how accurately 
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Figure 13: The sensitivity of the various Higgs pair production processes to the trilinear 
SM Higgs self-coupling at different cm. energies. The left panels display the total cross 
sections, the right panels display the ratio between the cross sections at a given k = 
^HHH / ^^HHH ^'^^ cross scctions at K = l. 



the respective Higgs pair production process has to be measured in order to extract the 
SM trihnear Higgs self-couphng with a certain accuracy. For example the gluon fusion 
cross section has to be measured with an accuracy of ~ 50% at = 8 TeV in order to 
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be able to extract the trilinear Higgs self-coupling with an accuracy of 50%, as can be 
inferred from Fig. [13] upper left. Similar dicussions can be found in Refs. [5|[T8]. 

4 Prospects at the LHC 

As shown in the previous section, inclusive Higgs boson pair production is dominated by 
gluon fusion at LHC energies. Other processes, such as weak boson fusion, qq' — )■ qq'HH, 
associated production with heavy gauge bosons, qq' — j- VHH (V = W, Z) or associated 
production with top quark pairs, gg/qq — >■ tiHH, yield cross sections which are factors of 
10 - 30 smaller than that for gg — j- HH. Since at the LHC Higgs boson pair production 
cross sections are small, we concentrate on the dominant gluon fusion process. In the 
following, we examine channels where one Higgs boson decays into a b quark pair and 
the other Higgs boson decays into either a photon pair, gg — )• HH — 6677, into a r 
pair, gg — )■ HH — )■ bbT~^T~, or into an off-shell W boson pair, gg HH — )• bbW*W*. 
Following the Higgs Cross Section Working Group recommendations |[32], we assume a 
branching ratio of 57.7% for a 125 GeV Higgs decaying into b quarks, 0.228% for the 
Higgs decaying into a photon pair, 6.12% for the Higgs decaying into a r pair and 21.50% 
for the Higgs decaying into off-shell W* bosons. 

At the time of the analysis, no generator existed for the signal process, but the matrix 
element for Higgs pair production in the gluon fusion channel was available in the fortran 
code HPAIR [18lll2]. In parallel to the approach used by the program described in [581159]. 
the HPAIR matrix element was added to Pythia 6 [60]. It has been verified that the cross 
sections produced by HPAIR and by the Pythia 6 implementation are consistent. 

All tree-level background processes are calculated using Madgraph 5 [61]. Signal and 
background cross sections are evaluated using the MSTW2008 parton distribution func- 
tions [SO] with the corresponding value of Ofg at the investigated order in perturbative 
QCD. The effects of QCD corrections are included in our calculation via multiplicative 
factors which are summarized in the following subsections and have been introduced in 
section 3. 

4.1 Kinematical distributions of gg — > HH 

In this subsection the characteristic distributions of the gluon fusion process gg — HH are 
studied for several observables. In Fig. [TH we show for each of the two final Higgs bosons 
the normalized distributions of the transverse momentum Pt,h and the pseudorapidity 
rjH, as well as the invariant mass Mhh, the helicity angle which is the angle between 
the off-shell Higgs boson, boosted back into the Higgs boson pair rest frame, and the Higgs 
boson pair direction, and the rapidity i/hh of the Higgs boson pair. The distributions of 
each observable are shown for different values of the trilinear Higgs coupling X/\sm = 
(green curve), 1 (red curve) and 2 (blue curve). We also include in the plots the typical 
background qq — t- ZH coming from the Higgs boson itself (black curve), the Z boson 
faking a Higgs boson. 

The Higgs bosons from inclusive Higgs pair production are typically boosted, as we 
can see in the upper left plot of Fig. [Ml where the Pt,h distributions reach their maximum 
for Pt,h ~ 150 GeV. For X/Xsm = 2, the triangle diagram interferes destructively with 
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[TeV] 


HH bb-fj 


ttH 


ZH 


14 


1.88 1.0 


1.10 


1.33 



Table 6: K -factors for gg -> HH, bb'yy, ttH and ZH production at y/s = 14 TeV JM/. 
The Higgs boson mass is assumed to be Mh = 125 GeV. 



the box diagram, which explains the dip in the Pt,h distribution. This high transverse 
momentum spectrum can also be interpreted in terms of the low pseudorapidity of the 
two Higgs bosons which have a typical symmetric distribution with the maximum around 
zero, see Fig. [14] upper right. The qq — ?■ ZH background has a completely different 
topology with less boosted Higgs and Z bosons, Pt,h/z ~ 50 GeV, with pseudorapidity 
of order \riH/z\ ~ 2 as can be seen in the upper left and right plots of Fig. [TH The middle 
left plot of Fig. [13] displays the distributions of the invariant mass of the Higgs boson 
pair. For the signal the typical value is M^h ^ 400 GeV to be compared to a lower 
value of MzH ^ 250 GeV for the ZH background. Also note that an important depletion 
appears in the signal for X/Xsm = 2 caused by the destructive interference between the 
triangle and box contribution. Due to the Higgs boson scalar nature a known discriminant 
observable is the angle O'^fj [62] . The middle right plot in Fig. [H] shows that signal and 
ZH background have similar distributions thus making this observable less discriminant 
than others described before but still efficient for some specific background, as we will see 
in the following. From the bottom plot of Fig. [141 it can be inferred that the rapidity 
distribution of the Higgs pair is narrower for the signal than for the ZH background. 

In the following the different decay channels HH — )■ 6677, HH — ?■ bbrf and HH — j- 
bbW^W~ will be investigated in more detail. 

4.2 The 6677 decay channel 

In this subsection, the 6677 final state for the production of two Higgs bosons with a mass 
of 125 GeV at ^/s = 14 TeV is investigated. Earlier studies can be found in Ref. |36] . 
The calculation of the signal, pp — ?■ HH — )■ 6677, is performed as described above by 
incorporating the matrix element extracted from the program HPAIR into Pythia 6. We 
include the effects of NLO QCD corrections on the signal by a multiplicative factor, 
Knlo = 1-88, corresponding to a 125 GeV Higgs boson and a cm. energy of 14 TeV. 
Here we set the factorization and renormalization scales equal to Mhh- The background 
processes, which have been generated, are the QCD process 6677 and the associated 
production of a Higgs boson with a ti pair, ttH , with the Higgs boson subsequently 
decaying into a photon pair and the top quarks decaying into a W boson and a b quark as 
well as single Higgs-strahlung ZH with the Higgs boson decaying into 77 and the Z boson 
decaying into hb. The QCD corrections have been included by a multiplicative i^-factor 
applied to the respective LO cross section. All i^-factors are taken at NLO except for 
the single Higgs-strahlung process which is taken at NNLO QCD and the case of 6677 in 
which no higher order corrections are taken into account. The various /'^-factors are given 
in Table [6] and taken from Ref. [33]. The factorization and renormalization scales have 
been set to Mhh for the signal and values specific for each process for the backgrounds. 
In this analysis, the signal and background processes are generated with exclusive 
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Figure 14: Normalized distributions of Pt^h, Vh, ^hh, (^hh ^'^^ Uhh for different values 
of the trilinear Higgs coupling in terms of the SM coupling, X/Xsm = 0, 1, 2. 

cuts. The basic acceptance cuts are motivated by the fact that the b quark pair and the 
photon pair reconstruct the Higgs mass according to the resolutions expected for ATLAS 
and CMS. 

In detail, we veto events with leptons having soft transverse momentum pT,e > 20 GeV 
and with a pseudorapidity \r]e\ < 2.4 in order to reduce the tiH background. Furthermore 
we also veto events with QCD jets prjet > 20 GeV and with a pseudorapidity \r]jet\ < 2.4 
to further reduce the tiH background. We require exactly one b quark pair and one 
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photon pair. The b quark pair is restricted to have pT,b > 30 GeV, |?7b| < 2.4 and 
AR{b,b) > 0.4, where AR{b,b) denotes the isolation of the two b quarks defined by 
the distance AR = (Ar/)^ + (A0)^ in the pseudorapidity and azimuthal angle plane 
(?7,0). We consider the fe-tagging efficiency to be 70%. The photon pair has to fulfill 
Pt,7 > 30 GeV, \r]^\ < 2.4 and AR{'j,j) > 0.4. The two reconstructed Higgs bosons, 
from the b quark pair and from the photon pair, have to reproduce the Higgs boson mass 
within a window of 25 GeV, 112.5 GeV < < 137.5 GeV, and a window of 10 GeV, 
120 GeV < < 130 GeV, respectively. We require additional isolations between the 
b quarks and the photons being Ai?(7, b) > 0.4. 

Based on the distributions shown in Fig. [TSl apart from the acceptance cuts we have 
applied more advanced cuts for this parton level analysis. We first require the recon- 
structed invariant mass of the Higgs pair to fulfill Mhh > 350 GeV. Furthermore we 
remove events which do not satisfy Pt.h > 100 GeV. We also constrain the pseudorapid- 
ity of the two reconstructed Higgs bosons, |?7h| < 2, and the isolation between the two b 
jets to be AR{b,b) < 2.5. 





AR(b,b) 

Figure 15: Normalized signal and backgrounds distributions for of Pt^h, M^h (md Rbf, in 
the bb^^ channel. 

The results are collected in Table [71 The local decrease of the sensitivity between the 
cut on Mhh and the cut on Pt,h is explained by the fact that we accept to have a reduced 
sensitivity locally during the chain of cuts in order to enhance the final significance. In the 
case described in this section a cut on Pt,h alone reduces the sentivity as does a cut on -qn 
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HH 6677 tt77 ZH S/B S/Vb 



Cross-section NLO [fb] 


8.92 X 10" 


-2 


5.05 X 10^ 


1.39 




3.33 X 10- 


-1 


1.77 X 10- 


-5 


6.87 X 10-2 


Reconstructed Higgs from bs 


4.37 X 10- 


-2 


4.01 X 10^ 


8.70 X 10- 


'2 


1.24 X 10- 


-3 


1.09 X 10- 


-4 


1.20 X 10-1 


Reconstructed Higgs from 7s 


O.UO X iU 


-2 


1. (O 


/1 w 1 n- 
z.4o X lU 


-2 


0. ( X iU 


-4 


i.oy X lu 


-2 


1 0/1 

1 .z4 


Cut on Mhh 


2.73 X 10- 


-2 


3.74 X 10-2 


7.45 X 10- 


-3 


1.28 X 10- 


-4 


6.07 X 10- 


-1 


7.05 


Cut on Ptm 


2.33 X 10- 


-2 


3.74 X 10-2 


5.33 X 10- 


-3 


1.18 X 10- 


-4 


5.44 X 10- 


1 


6.17 


Cut on TjH 


2.04 X 10- 


-2 


1.87 X 10-2 


3.72 X 10- 


3 


9.02 X 10- 


-5 


9.06 X 10- 


-1 


7.45 


Cut on Ai?(6, h) 


1.71 X 10- 


-2 


0.00 


3.21 X 10- 


3 


7.44 X 10- 


^5 


5.21 




16.34 


"Detector level" 


1.56 X 10- 


-2 


0.00 


8.75 X 10- 


-3 


8.74 X 10- 


-3 


8.92 X 10- 


1 


6.46 



Table 7: Cross-section values of the HH signal and the various backgrounds expected at 
the LHC at y/s = 14 TeV, the signal to background ratio S/B and the significance S/\fB 
for j C = 3000 /^-^ in the 6677 channel after applying the cuts discussed in the text. 



alone, but the first cut actually improves the discrimination between the signal and the 
background in the pseudorapidity distribution, hence allowing for a larger improvement 
when applying the rju cut just after the Pt,h cut. Eventually all the cuts allow for an 
improvement of the significance by two orders of magnitude, that is the ratio of signal 
events S over the square root of background events B, S/ \fB. The final value for S/y/B 
is 16.3 for an integrated luminosity of j£ = 3000 fb~^, corresponding to 51 signal events. 
Therefore this channel seems promising. 

A realistic assessment of the prospects for measuring the signal in the 6677 final 
state depends mostly on a realistic simulation of the diphoton fake rate due to multijet 
production, which is the dominant background in such an analysis. Our first parton level 
study gives a rough idea of how promising the 6677 channel is. 

In the following we perform a full analysis including fragmentation and hadronisation 
effects, initial and final state radiations by using Pythia 6 . 4 in order to assess more 
reliably whether the promising features of this channel survive in a real experimental 
environment. All the events are processed through Delphes [63], the tool which is used 
for detector simulation. For the jet reconstruction we use the anti-fcr algorithm with 
a radius parameter of i? = 0.5. We still consider a 6-tagging efficiency of 70%. We 
keep the same acceptance cuts as before except for the transverse momentum of the 
reconstructed b jet and photon which we increase up to Pt,6/7 > ^0 GeV. We also enlarge 
the window for the reconstructed Higgs boson coming from the b quark pair, by requiring 
100 GeV < Mf^i < 135 GeV. We select events with exactly two reconstructed b jets and 
two photons. 

The final result is displayed in the last line of Table [71 The final significance S/ y/B for 
this simulation has decreased to 6.5 for an integrated luminosity of 3000 fb~^, correspond- 
ing to 47 events. Though low, the significance nevertheless is promising enough to trigger 
a real experimental analysis as can be performed only by the experimental collaborations 
and which is beyond the scope of this work. 
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[TeV] 


HH hbr+T" 


tt 


ZH 


14 


1.88 1.21 


1.35 


1.33 



Table 8: K-factors for gg HH, bbrr {6^, tt J6^ and ZH production JM^ at ^ = 
14 TeV. The Higgs boson mass is assumed to be Mh = 125 GeV. 



4.3 The hhrr decay channel 

The hbrf decay channel is promising for low mass Higgs boson pair production at the 
LHC and has been previously studied in Ref. [37] . An important part of this analysis 
will depend on the ability to reconstruct the b quark pair and the r pair. As the real 
experimental assessment of such a reconstruction is beyond the scope of our work we will 
perform in the following a parton level analysis, assuming a perfect r reconstruction. The 
analysis thus respresents an optimistic estimate of what can be done at best to extract 
the trilinear Higgs self-coupling in the hbrf channel. 

We consider the two QCD-QED continuum background final states hbrf and bbrv^-rVr 
calculated at tree-level. The bbrUrTUr final state background dominantly stems from ti 
production with the subsequent top quark decays t — )■ bW and W fUr. We also 
include backgrounds coming from single Higgs production in association with a Z boson 
and the subsequent decays H ^ rf and Z ^ bb or H ^ bb and Z ^ rf. The effects of 
QCD corrections are included in our calculation via multiplicative JiT-factors which are 
summarized in Table All i^-factors are taken at NLO except for the single Higgs- 
strahlung process which is taken at NNLO QCD. The factorization and renormalization 
scales have been taken at Mhh for the signal and at values specific for each background 
process. 

Concerning the choice of our cuts, we demand exactly one b quark pair and one r 
pair. The b quark pair is required to fulfill pT,b > 30 GeV and |?7fe| < 2.4. We assume 
the 6-tagging efficiency to be 70% and the r-tagging efficiency to be 50% and we neglect 
fake rates in both cases. The r pair has to fulfill pt,t > 30 GeV and |?7t-| < 2.4. The 
reconstructed Higgs boson from the b quark pair is required to reproduce the Higgs mass 
within a window of 25 GeV, 112.5 GeV < M^j < 137.5 GeV. The reconstructed Higgs 
boson from the r pair needs to reproduce the Higgs mass within a window of 50 GeV, 
100 GeV < Mrf < 150 GeV or within a window of 25 GeV, 112.5 GeV < M^f < 
137.5 GeV, in a more optimistic scenario. In addition to these acceptance cuts we also 
add more advanced cuts for this parton level analysis, based on the distributions shown in 
Fig. [16] and in a similar way as what has been done in the previous 6677 analysis. We first 
demand the invariant mass of the reconstructed Higgs pair to fulfill M^h > 350 GeV. In 
addition, we remove events which do not satisfy Pt^h > 100 GeV. We do not use a cut on 
the pseudorapidity of the reconstructed Higgs bosons in this analysis as it would reduce 
the significance. 

The different results of our parton level analysis are summarized in Table |9l The cuts 
allow for an improvement of two orders of magnitude in the significance S/ a/B. The final 
S/ \/B is 6.71 for an integrated luminosity oi JC = 3000 fb~^, corresponding to 330 signal 
events. We then conclude that this channel is promising. In the last line we reproduce 
our result for the optimistic requirement of 112.5 GeV < Mrr < 137.5 GeV leading to 
the final significance S/\fB = 9.36 for an integrated luminosity of 3000 fb~^. Already 
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Figure 16: Normalized distributions of Pt,h, Mhh O'^d tjh for signal and backgrounds in 
the bbrf channel. 

for a planned mid-term integrated luminosity of 300 fb~^ at 14 TeV the expectations 
are promising with 33 signal events and a significance S/\fB = 2.96 in the optimistic 
scenario. 

4.4 The bbW^W^ decay channel 

The analysis in this channel is difficult as the leptonic W boson decays lead to missing 
energy in the final state. Consequently, one of the two Higgs bosons cannot a priori 
be reconstructed equally well as the other Higgs boson, thus reducing our capability to 
efficiently remove the background with the canonical acceptance cuts previously applied 
in the other decay channels. This channel with one lepton plus jets final state has been 
studied in jnZlEH]- We only consider here the decay W — )■ iui {i = e, //) with a branching 
ratio of 10.8%. 

We take into account the continuum background which contains all processes with 
the bbiueiui final states at tree-level, for example qq/gg — )■ b*b* — > 'jbZb — bbiiZ with 
the subsequent splitting Z — )■ lygiyeii- We proceed in a similar manner as in the previous 
analyses. We generate the signal and the backgrounds with the following parton-level 
cuts. We require that the b quarks fulfill pT,b > 30 GeV and \r]b\ < 2.4. We consider the 
6-tagging efficiency to be 70%. The leptons have to fulfill pt/ > 20 GeV and \r]i\ < 2.4. 
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HH 




bbr^T 


bbr'^T Vt^Vt- 


ZH 




S/B 




s/Vb 


Cross-section NLO [fbl 

L J 


2.47 




2.99 X 10'' 


8.17 X 10^ 


2.46 X 10 


1 


6.48 X 10- 


-5 


6.93 X 10-1 


PTonstriif'tpr] TTie^'P's from ts 


2.09 X lO'" 


-1 


8 35 X 10^ 


1.58 X 10^ 


5 70 X 10- 


-1 


8 63 X 10" 


-4 


7 36 X IQ-i 


Reconstructed Higgs from 6s 


1.46 X 10" 


-1 


6.34 X 10-1 


1.43 X 10^ 


3.75 X 10- 


-2 


9.75 X 10- 


-3 


2.07 


Cut on Mhh 


1.30 X 10- 


-1 


1.37 X 10-1 


1.74 


1.26 X 10- 


-2 


6.88 X 10" 


'2 


5.18 


Cut on Pt,h 


1.10 X 10" 


-1 


7.80 X 10-2 


7.17 X 10-1 


1.15 X 10- 


-2 


1.36 X 10- 


1 


6.71 


With 112.5 < Mrr < 137.5 


1.10 X 10" 


-1 


3.41 X 10-2 


3.76 X 10-1 


3.15 X 10- 


-3 


2.67 X 10- 


-1 


9.37 



Table 9: Cross-section values of the of HH signal and the various backgrounds expected at 
the LHC at y/s = 14 Te V, the signal to background ratio S/ B and the significance S/ \fB 
for j jC = 3000 fb^^ in the bbrf channel after applying the cuts discussed in the text. 



The reconstructed Higgs boson from the b quark pair has to reproduce the Higgs boson 
mass within a window of 25 GeV, 112.5 GeV < M^^ < 137.5 GeV. We also require that 
the missing transverse energy respects E^^'""^ > 20 GeV. 

As done in the previous subsections, we also add more advanced cuts for this parton 
level analysis, based on the distributions shown in Fig. [T71 The distributions on the 
upper left of Fig. [17] correspond to the transverse mass of the lepton pair, being defined 
as Mt = ^2p^rfE^'''{l - cosA0(E^^^^££)), where A(j){EJp^'' , M) is the angle between 
the missing transverse momentum and the transverse momentum of the dilepton system. 
The distribution of the signal has an endpoint at the value of Mh- The distributions 
on the upper right of Fig{T7] represent the angle between the two leptons projected on 
the transverse plane, Ac^^^^j. The angle is reduced for the signal compared to the broad 
distribution of the background. The last distributions on the bottom of Fig. [T7] display 
the projected missing transverse energy E^^^'^ = i?™"** sin A0(_E™'**, £) for A(/)(£'^*^*, i) < 
7r/2, where A(^(£'™^*,£) is the angle between the missing transverse momentum and 
the transverse momentum of the nearest lepton candidate. If A0(i?™^*,£) > 7r/2, then 
j^miss _ j^mtss^ The signal distribution is shifted to the left compared to the background 
distribution. 

We first require the transverse mass of the lepton pair to be Mt < 125 GeV. We then 
remove events which do not satisfy A^^^^j < 1.2 and we also add a constraint on the angle 
between the two leptons, A^^^^g < 1.0. We demand the missing transverse energy to fulfill 
^jniss ^ GeV and the projected energy to satisfy EJp^'' < 80 GeV. Note that the 
j^miss distribution displayed in Fig. [T7] is obtained after the acceptance cuts having been 
applied but before the advanced cuts. The cuts on Mr, ^(pi^i^, ^Ghi2 ^^'^ modify 
this distribution and explain why the i?™'^^ cut, which would seem not to be efficient, 
actually improves the significance. 

The results for the LHC at y/s = 14 TeV are summarized in Table [TOl While the cuts 
allow for an improvement of the significance S/ \/B by about one order of magnitude, we 
are still left with a a very small signal to background ratio. Thus, this channel using the 
final states considered here is not very promising. 
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Figure 17: Normalized distributions of Mt, ^<Phi2 and projected missing transverse energy 
T^miss signal and background channels in the bblih'ij2^h i^^^^ states of the bbW^W~ 
channel. 





HH 




bbliUij2J^l^ 


S/B 








Cross-section NLO [fb] 


3.92 X 10-^ 


-1 


2.41 X 10^ 


1.63 X 10' 


-5 


1.38 X 10" 


-1 


Reconstructed Higgs from bs 


6.18 X 10" 


-2 


1.89 X 10^ 


3.27 X 10" 


-4 


2.46 X 10" 


-1 


Cut on Mt 


6.18. X 10" 


-2 


1.19 X 10^ 


5.19 X 10" 


-4 


3.10 X 10" 


-1 


Cut on A0£^£2 


5.37 X 10' 


-2 


6.96 X 10^ 


7.72 X 10" 


-4 


3.53 X 10" 


-1 


Cut on A9i-^£^ 


5.17 X 10" 


-2 


5.65 X IQi 


9.15 X 10" 


-4 


3.77 X 10" 


-1 


Cut on EJp^'' 


8.41 X 10" 


"3 


3.77 X 10-1 


2.22 X 10" 


-2 


7.50 X 10" 


-1 


Cut on EJp^""^ 


4.59 X 10" 


-3 


2.70 X 10-2 


1.70 X 10" 


-1 


1.53 





Table 10: Cross-section values of the HH signal and the considered background expected 
at the LHC at ^/s = 14 TeV, the signal to background ratio S/B and the significance 
S/\fB for JC=3000 fb^^ in the bbW^W~ channel after applying the cuts discussed in 
the text. 
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5 Conclusions 



In this paper we have disussed in detail the main Higgs pair production processes at the 
LHC, gluon fusion, vector boson fusion, double Higgs-strahlung and associated production 
with a top quark pair. They allow for the determination of the trilinear Higgs self-coupling 
^HHHi which represents a first important step towards the reconstruction of the Higgs 
potential and thus the final verification of the Higgs mechanism as the origin of electroweak 
symmetry breaking. We have included the important QCD corrections at NLO to gluon 
fusion and vector boson fusion and calculated for the first time the NNLO corrections to 
double Higgs-strahlung. It turns out that the gluon initiated process to ZHH production 
which contributes at NNLO is sizeable in contrast to the single Higgs-strahlung case. We 
have discussed in detail the various uncertainties of the different processes and provided 
numbers for the cross sections and the total uncertainties at four cm. energies, i.e. 8, 14, 
33 and 100 TeV. It turns out that they are of the order of 40% in the gluon fusion channel 
while they are much more limited in the vector bosons fusion and double Higgs-strahlung 
processes, i.e. below 10%. Within the SM we also studied the sensitivities of the double 
Higgs production processes to the trilinear Higgs self-coupling in order to get an estimate 
of how accurately the cross sections have to be measured in order to extract the Higgs 
self-interaction with sufficient accuracy. 

In the second part of our work we have performed a parton level analysis for the domi- 
nant Higgs pair production process through gluon fusion in different final states which are 
6677, bbrf and bbW^W~ with the W bosons decaying leptonically. Due to the smallness 
of the signal and the large QCD backgrounds the analysis is challenging. The hhW^W~ 
final state leads to a very small signal to background ratio after applying acceptance and 
selection cuts so that it is not promising. On the other hand, the significances obtained 
in the 6677 and hhrf final states after cuts are ~ 16 and ~ 9, respectively, with not 
too small event numbers. They are thus promising enough to start a real experimental 
analysis taking into account detector and hadronization effects, which is beyond the scope 
of our work. Performing a first simulation on the detector level for the 6677 state shows, 
however, that the prospects are good in case of high luminosities. 
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